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The absolute branching fraction for the decay D + — > fi + u has been directly measured based on 
a data sample of about 33 pb _1 collected around yfs = 3.773 GeV with the BEST! detector at the 
BEPC collider. A total of 5321 ± 149 ± 160 D~ mesons are reconstructed in nine hadronic decay 
modes. In the system recoiling against these singly tagged D~ mesons, 2.67 ± 1.74 purely leptonic 
decay events of D + — > pi + f M are observed. Those yield the branching fraction of BF(D + — > 
^i + f M ) = (0.122l° Q53 ± 0.010)%, and a corresponding value of the pseudoscalar decay constant 
f D+ = (371±i 20 ± 25) MeV. 
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son. The virtual W + boson is produced in the anni- 
hilation of the c and d quarks. The decay rate of this 
process is determined by the wavefunction overlap of the 
two quarks at the origin, and is parametrized by the D + 
decay constant, fu+. Fig. 1 shows the decay diagram 
for the Cabibbo-suppressed purely leptonic decay of the 
D + meson. The decay width of the D + — > l + vi is given 
by the formula 



BES Collaboration 141] measured f D + = (3001^112) 

MeV based on one event of D + — ► from the data 

collected at 4.03 GeV with the BES-I detector at the 
BEPC collider. 

In this Letter, we report a direct measurement of the 

branching fraction for the decay D + — > and deter- 
mination of the decay constant f D + . 



T(D+ - l+u,) = ^± | Vcd | 2 mfm D+ (l - -4~) , 
8vr V m- D+ J 

(1) 

where Gf is the Fermi coupling constant, V c d is the c — > <i 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element 0, 
to/ is the mass of the lepton, and m D + is the mass of the 
D + meson. 




FIG. 1: The decay diagram for D — > I V\. 

fi)+ is an important parameter. However, it is diffi- 
cult to measure f D + due to the fact that the D + — ► l + vi 
is a Cabibbo-suppressed decay process. There are some 
theoretical calculations to estimate the value of fn+ ■ Al- 
though predictions for f D + vary significantly from 90 to 
360 MeV (3, the predictions for the ratios fu+ : f D + : fg 
are more reliable, where f D + and fg are the decay con- 
stants for the Dg and the B mesons, respectively. The 
meson decay constants play an important role in extract- 
ing some interesting physics quantities from diverse mea- 
surements. For example, fg relates the measurement of 
the BB mixing p| ratio to CKM matrix elements. At 
present it is not possible to determine fg experimentally 
from the purely leptonic B decay, so theoretical calcu- 
lations of fg must be used. Hence, the calculations of 
fs are of considerable importance. With the predictions 
for the ratios fjj+ : f D + : fg, measurements of fu+ and 
f D + provide checks on some theoretical calculations of 
the decay constants and help discriminate among differ- 
ent models and improve the reliability of estimates of fg . 

To date, there are eight experimental measurements of 
f D + from the WA75 j|, CLEO % BES 0, E653 0, 
L3 S @, BEATRICE pjj, OPAL Q3 and ALEPH [12 
groups. For D+ -> fi+v^, the MARK-III Collabora- 
tion 0| set a branching fraction upper limit of 0.07% 
fcorresDondine to f n + < 290 MeV at 90% C.L.I. The 



II. THE BES-II DETECTOR 

The BES-II is a conventional cylindrical magnetic de- 
tector that is described in detail in Ref. [15j. A 12-layer 
vertex chamber (VC) surrounding the beryllium beam 
pipe provides input to the event trigger, as well as co- 
ordinate information. A forty-layer main drift chamber 
(MDC) located just outside the VC yields precise mea- 
surements of charged particle trajectories with a solid an- 
gle coverage of 85% of An; it also provides ionization en- 
ergy loss (dE/dx) measurements which are used for parti- 
cle identification. Momentum resolution of 1.7% y/l + p 2 
(p in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering electrons are obtained for the data taken at 
= 3.773 GeV. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) 
of charged particles with a resolution of about 180 ps 
for electrons. Outside the TOF, a 12 radiation length, 
lead-gas barrel shower counter (BSC), operating in self- 
quenching streamer mode, measures the energies of elec- 
trons and photons over 80% of the total solid angle with 
an energy resolution of cte/E = 0.22/v^E (E in GeV) 
and spatial resolutions of cr^ = 7.9 mrad and <jz = 2.3 
cm for electrons. A solenoidal magnet outside the BSC 
provides a 0.4 T magnetic field in the central tracking re- 
gion of the detector. Three double-layer muon counters 
instrument the magnet flux return, and serve to iden- 
tify muons of momentum greater than 500 MeV/c. They 
cover 68% of the total solid angle. 



III. DATA ANALYSIS 

The data used in the analysis were collected with the 
BES-II detector at the BEPC collider. A total integrated 
luminosity of about 33 pb -1 was taken at and around the 
center-of-mass energy of 3.773 GeV. Those are just above 
the threshold of e + e _ — > DD and below the threshold 
of e + e~ — > DD . Thus, if a D~ meson decay is fully 
reconstructed (This is called a singly tagged D~ meson), 
a D + meson must exist in the system recoiling against 
the singly tagged D~ meson. From the singly tagged D~ 
event sample, the events of the decay D + — > [i^v^ can 
be well selected in the recoiling system. Therefore, the 
absolute branching fraction for the decay D + — > p^v^ 
can be well measured and the decay constant f D + can be 
determined. 
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A. Singly tagged D event sample 

1. Events selection 



The D~ meson is reconstructed in the nine hadronic 
decay modes of K + ir~ir~, K°ir~, K°K~, K + K~tt~ , 
K°ir~ir~ir + , A" 7r~7T , K + ir~ir~ir°, A" + 7r + 7r~7r~7r~ and 
7r + 7r~7r~. Events which contain at least three recon- 
structed charged tracks with good helix fits are selected. 
In order to ensure the well-measured 3-momentum vec- 
tors and the reliably charged particle identification, the 
charged tracks used in the single tag analysis are required 
to be within \cos8\ <0.85, where is the polar angle of 
the charged track. All tracks, save those from Kg de- 
cays, must originate from the interaction region, which 
require that the closest approach of the charged track in 
the xy plane is less than 2.0 cm and the absolute z posi- 
tion of the track is less than 20.0 cm. Pions and kaons are 
identified by means of TOF and dE/dx measurements. 
Pion identification requires a consistency with the pion 
hypothesis at a confidence level (CL n ) greater than 0.1%. 
In order to reduce misidentification, a kaon candidate is 
required to have a larger confidence level (CLx) for a 
kaon hypothesis than that for a pion hypothesis. The ir° 
is reconstructed in the decay of ir° — > 77. To select good 
photons from the decay of 7r°, the energy of a photon 
deposited in the BSC is required to be greater than 0.07 
GeV 0|, and the electromagnetic shower is required to 
start in the first 5 readout layers. In order to reduce back- 
grounds, the angle between the photon and the nearest 
charged track is required to be greater than 22° and 
the angle between the direction of the cluster develop- 
ment and the direction of the photon emission to be less 
than 37° O. 

For the single tag modes of D~ — > A" + 7r + 7r~7r~7r _ and 
D~ — > 7r + 7r~7r~ , backgrounds are further reduced by re- 
quiring the difference between the measured energy of 
the D~ candidate and the beam energy to be less than 
70 and 60 MeV, respectively. In addition, the cosine of 
the D~ production angle relative to the beam direction 
is required to be \cos8 D - \ < 0.8. 



2. Single tag analysis 

For each event, there may be several different charged 
track (or charged and neutral track) combinations for 
each of the nine single tag modes. Each combination 
is subject to a one-constraint (1C) kinematic fit requir- 
ing overall event energy conservation and that the un- 
measured recoil system has the same invariant mass as 
the track combinations. Candidates with a fit probabil- 
ity P(x 2 ) greater than 0.1% are retained. If more than 
one combination satisfies P(x 2 ) > 0-1%, the combination 
with the largest fit probability is retained. For the single 
tag modes with a neutral kaon and/or neutral pion, one 
additional constraint kinematic fit for the Kg — > ir + ir~ 
a/nd/or 7r° — > •yy hvnothesis is nerformed. senaratelv. 



The resulting distributions in the fitted invariant 
masses of mKnn (m — or 1 or 2 and n = 1 or 2 or 3 or 4) 
combinations, which are calculated using the fitted mo- 
mentum vectors from the kinematic fit, are shown in 
Fig. |21 The signals for the singly tagged D~ mesons 
are clearly observed in the fitted mass spectra. A max- 
imum likelihood fit to the mass spectrum with a Gaus- 
sian function for the D~ signal and a special background 
function^ to describe backgrounds yields the number 
of the singly tagged D~ events for each of the nine modes 
and the total number of 5321 ± 149 ± 160 reconstructed 
D~ mesons, where the first error is statistical and the 
second systematic obtained by varying the parameteri- 
zation of the background. The curves of Fig. [2] give the 
best fits to the invariant mass spectra. In the fits to the 
mass spectra, the standard deviations of the Gaussian 
signal functions for Fig. Otg) and Fig. |2Ih) are fixed at 
4.27 MeV and 2.16 MeV, respectively. These standard 
deviations are obtained from Monte Carlo sample. All 
other parameters are left free in the fit. 
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FIG. 2: The distributions of the fitted invariant 
masses of (a) K^~tt~ tt~ , (b) K°tv~ , (c) K°K~ , (d) 
K+K--K-, (e) K°-k-tv-tt + , (f) A°7r-7r°, (g) K + ir~ ir~ ir° , 
(h) K + 7r + tv~ ir~ tv~ , and (i) tv~tv~tt + combinations; (j) the 
9 modes combined together; the curves are the best fits de- 
scribed in the text. 
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B. Events of IT 



1. Muon Identification 



To identify the muon from the decay D + — > fj vv^, the 
charged track in the recoil system of the D~ tag in each 
of the events as shown in Fig. [21 is examined for muon 
identification requirements. The muon identification re- 
quires: 

1. The charged track satisfies | cos6>| < 0.68, where 9 
is the polar angle of the charged track. 

2. There must be hits in the muon system and the hits 
must well associate with the charged track extrap- 
olated from the track reconstructed in the MDC. 
The required number of the hits in the muon sys- 
tem is momentum dependent. 

3. The muon candidate with the transverse momen- 
tum of greater than 0.7 GeV/c is required to hit at 
least two layers of the muon system. 

A charged track satisfying the 3 requirements is identified 
as a muon. However, a small fraction of the pions which 
can punch through the muon system could be misidenti- 
fied as muons. 

To estimate the average probability of misidentify- 
ing a pion as a muon with momentum between 0.785 
and 1.135 GeV/c, we select the pions from the decay 
J '/ip — > lott + tt~ in the data taken at the center-of-mass 
energy of 3.097 GeV with the BES-II detector. The pion 
which satisfies |cos#,r| < 0.68 (6^ is the polar angle of 
the pion) and with momentum in the region from 0.785 
to 1.135 GeV/c is checked for whether it satisfies the 
muon selection requirements. A total of 10657 pions 
from the J/ip — ► ujtt + tt~ events satisfy the pion selec- 
tion criteria and 259 of them are misidentified as muons. 
Those give the averaged misidentification probability of 
0.024 ± 0.002. 



2. Candidates for D + 



Candidate events for the decay D + — ► ^ 
lected from the surviving charged track in the system 
recoiling against the singly tagged D~ mesons. To select 
the D + — > /i + t^, it is required that there be a single 
charged track originating from the interaction region in 
the recoil system of the D~ tag and the charged track is 
identified as a muon with charge opposite to the charge 
of the tagged D~. For the candidate event, no extra 
good photon which is not used in the reconstruction of 
the singly tagged D~ meson is allowed to be present in 
the event. Since there is a missing neutrino in the purely 
leptonic decay event, the event should be characteristic 
with missing energy E m i SS and missing momentum P m i SS 
which are carried by the neutrino. For the purely leptonic 
decav event, the reconstructed U~,t.. which is defined as 



the difference between the E m i SS and the P m iss should 
be close to zero. Fig. OJa) shows the distribution of the 
Umiss for the Monte Carlo events of e + e~ — > D + D~, 



where D 



K + tt 7r and D + 



shows the same distribution of the U„ 



iu- Fig. EIb) 

for the Monte 



Carlo events of D~ 



A' H 



tt u and £> H 



The distribution of the reconstructed muon momentum 
selected from the Monte Carlo events of D + — > 
is shown in Fig. 0] where the interval between the two 
dashed lines is defined as the selected muon momentum 
region for the events of D + — > /! + z^. 
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FIG. 3: The distributions of the U m iss calculated for the 

Monte Carlo events of (a) D + — + versus the tags of 

D~ — > K + -k~it~ and (b) D + — » versus the tags of 

D~ -> A' + TT-TT-TT . 
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FIG. 4: The distribution of the reconstructed muon momen- 
tum selected from the Monte Carlo events of D + — > 
the interval between the two dashed lines shows the region of 
the selected muon momentum in the data analysis. 
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To select the purely leptonic decay events from the 
singly tagged D~ event sample, it is required that the 
Umiss of the candidate events should be within the 
±3cr[/ miss , region for the single tag mode(z), where 



is the standard deviation of the U„ 



distri- 



bution obtained from the Monte Carlo simulation for the 
event of D + — > versus the single tag mode(i) (i = 1 

is for K + tt~it~; i = 2 is for K°ir~... and i = 9 is 
for 7r + 7r~7r~ mode). A further criterion requires that 
the candidate muon momentum should be in the region 
from 0.785 to 1.135 GeV/c as shown in Fig. H Fig. 
shows the scatter-plot of the momentum of the candi- 
date muon versus the U m i SS , where the dots and star are 
for the tag modes of K + ir~ir~ and K + ir~ir~ir° , respec- 
tively. The solid (dashed) vertical lines give the ±3(ru mU , 
interval for the single tag mode of D~ — > K + ir~ir~ 
(D~ — > K + ir~ir~Tr°), while the two horizontal lines give 
the selected momentum region for the /i + from the purely 
leptonic decays of the D + mesons. The regions sur- 
rounded by the lines are defined as the selected signal 
regions for the two single tag modes. There are 3 events 
within the signal regions for the two single tag modes. 



TABLE I: Three candidates for D + purely leptonic decay. 

Event 1 2 3 

Tagging mode K + iv~ir~ K + ix~ iv~ ir° K + ir~ir~ 

Fitted mass [MeV/c 2 ] 1870.8 1876.9 1871.4 

Number of fi layer hits 2 2 2 

H + momentum [GeV/c] 0.974 0.981 0.919 

U mi3S [GeV] -0.093 -0.023 0.117 
Calculated momentum 

of neutrino [GeV/c] 1.000 1.007 0.843 



mass spectrum. Table [I] summarizes the characteristics 
of the 3 events. In each case, the measured mass of D~ 
meson, the U m i SS of the event and the momentum of the 
fi + agree well with the expected values. 
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FIG. 5: The scatter-plot of the muon momentum versus 
Umiss recoiling against the muon and the single tags (mKmn 
combinations) for the surviving purely leptonic decay candi- 
dates, where the dots (•) and the star (*) are for the single 
tag modes of K + tv~ty~ and K + ty~ tv~ it , respectively, those 
are from the data; the open squares and the open triangles 
are for the Monte Carlo background events from the Monte 
Carlo sample which is 27 times larger than the data, see text. 

Fig- EI a ) shows the distribution of the fitted masses of 
the mKnw combinations for the events which satisfy the 
selection criteria. The fitted masses of the 3 candidate 
events arc all in the well measured D~ sianal reeion in the 
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(b) 
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FIG. 6: The distributions of the mKmv (m — or 1 or 2 and 
n = 1 or 2 or 3 or 4) combinations for the events in which (a) 
the candidate events for D + — > fi + u^ are found in the system 
recoiling against the tags (mKrvn combinations) and (b) only 
one charged track is found in the system recoiling against the 
tags and the events satisfy the kinematic requirements for 
selecting the purely leptonic decay events, but these events 
do not satisfy muon identification requirements. 



3. Background subtraction 

Some non-purely leptonic events from the D + decays 
may also satisfy the selection criteria and are the back- 
around to the Durelv leDtonic decav events. These back- 
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ground events must be subtracted. The number of the 
background events is estimated by analyzing the Monte 
Carlo sample which is 27 times larger than the data. The 
Monte Carlo events are generated as e + e~ — * DD and 
the D and D mesons are set to decay to all possible final 
states according to the decay modes and the branching 
fractions quoted from PDG 01 except the purely leptonic 
decay mode under study. In Fig. [5] the open squares (□) 
and the open triangles (A) are for the background events 
from the Monte Carlo sample, where the □ and the A 
indicate that the background events are within and out- 
side of the ±3<Ju miss i interval, respectively. There are 9 
background events satisfying the selection criteria in the 
signal regions. This number of the background events is 
then normalized to the corresponding data set. A total 
of 0.33 ± 0.11 background events are estimated in the 3 
candidates for D + — ► n + v^. After subtracting this num- 
ber of the background events, 2.67 ± 1.74 signal events 
for D + — > decay are retained. 

The number of the background events can also be es- 
timated from the number of the singly tagged D~ events 
in which only one charged track in the recoiling system is 
found and the charged track is not identified as a muon. 
Fig.HJb) shows the distribution of the fitted masses of the 
mKniT combinations from the events. These events sat- 
isfy the kinematic requirements for selecting the purely 
leptonic decay events, but the charged track in the re- 
coiling system do not satisfy the muon identification re- 
quirements. There are 13 events in the D~ signal re- 
gions as shown in the hatched histogram. The average 
probability of misidentifying a pion as a muon discussed 
previously yields the number of the background events in 
the 3 purely leptonic decay candidates to be 0.31 ± 0.09, 
which is consistent with 0.33 ±0.11 estimated from the 
Monte Carlo sample 20] . 



IV. RESULT 

A. Monte Carlo Efficiency 

The efficiency for reconstruction of the purely leptonic 
decay D + — > after tagging the D~ meson is esti- 

mated by Monte Carlo simulation. A Monte Carlo study 
gives that the efficiency is e M +^ = (41.7 ± 1.1)%. 



N tag be the number of the singly tagged D~ mesons 
(Ntag = 5321 ± 149 ± 160), « M +„ M and n tag be the cor- 
responding expected numbers of the events, respectively. 
The likelihood function is then given by 

L = P{n^^,N^^)G(n tag ,N tag ), (2) 

with abbreviation 

r>t at \ K n ii+Vp) " -n + 



and 



G(n ta g, N tag ) 



-I t*f,» "t,,,,i 

e "top 



(3) 



where the a Ntag is the systematic uncertainty in the num- 
ber of the singly tagged D~ mesons ((TN tag = 160). 

For the observed numbers of the purely leptonic decay 
events and the singly tagged D~ mesons, the expected 
number of the purely leptonic decay events is given by 



— ntagBFe^+v^ + n b , 



(4) 



where BF is the purely leptonic decay branching fraction 
and rib is the number of the background events. The value 
of the likelihood function is obtained by integrating over 



L(BF) 



J L(BF,N^ +Uft )dN^^, 



(5) 



which is shown as a function of BF in Fig. [7| The 
maximum likelihood occurs at BF = 0.122%. Integrat- 
ing the function from the maximum position to — la 
(+lc) value corresponding to 68.3% of the total area 
below (above) the peak position yields the statistical er- 
ror to be -0.053% (+0.111%). The relative systematic 
uncertainty arises mainly from the uncertainties in the 
/i + identification (±5.0%), tracking efficiency (±2.0%), 
background subtraction (±5.6%) and U m i SS cut (±1.0%). 
Adding these uncertainties in quadrature gives the total 
relative systematic uncertainty to be ±7.8%. Finally, we 
obtain 



BF(D+ /*+!/„) 



(0.122^ ±0.010)%. 



B. Branching fraction 



To determine the purely leptonic decay branching frac- 
tion and the decay constant based on the observed num- 
bers of the purely leptonic decay events and the singly 
tagged D~ mesons, we build a likelihood function L, 
which is the product of the Poissonian probability func- 
tion for observation of the purely leptonic decay events 
and the Gaussian function for the number of the singly 
tagged D~ mesons. Let N^+y be the number of the 
observed Durelv leDtonic decav events (N..+,, = 3.01. 



C. Decay constant f D + 

The decay constant f D + can be obtained by inserting 
the mass of the muon, the mass of the D + meson, the 
CKM matrix element \V c d\, the Fermi coupling constant 
Gf and the lifetime of the D + meson into equation 
(1). Bv substituting BF(D + — > u + u,,) in terms of fn+ 
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0.002 0.004 0.006 

BF(D+ -» 

FIG. 7: The dependence of the normalized likelihood on the 
branching fraction for the decay of D + — > 

into equation (5), we obtain the relation of the likelihood 
function with fjj+ , 

L(f D +) = J Hfn^N^vJdN^^. (6) 

The most probable value of fo+ and its statistical error 
can be obtained by integrating over N IJj + v . The depen- 
dence of the likelihood function on the decay constant 
fu+ is shown in Fig. |SJ Following the procedure for ex- 
tracting the BF(D + — > we finally obtain 

fn+ = (37ltnl ± 25) MeV, 

where the first error is statistical and the second system- 
atic which arises mainly from the uncertainties in the 
measured branching fraction (±4.1%), the CKM matrix 
element \V c d\ (±5.4%), and the lifetime of the D + meson 
(±0.3%) @. These uncertainties are added in quadrature 
to obtain the total systematic error, which is ±6.8%. 

V. SUMMARY 

From the 5321 ± 149 ± 160 singly tagged D~ mesons, 
2.67 ± 1.74 purely leptonic decay events of D + — > 



are observed in the system recoiling against the singly 
tagged D~ mesons, resulting in the branching fraction 
of BF(D+ -> fi+Vf,) = (0.122±g;J|| ± 0.010)% and the 
decay constant of f D + = (371 1 \H ± 25) MeV. The mea- 
sured values are independent of the D + D~ cross section 
and luminosity. They do not require model-dependent as- 
sumptions. Thus, they are absolute measurements. The 
central value of fu+ is consistent with that measured us- 
ing the BES-I detector The measured value of f D + 
is also consistent with most theoretical predictions, which 
are in the range from 90 to 360 MeV [3j. 




/d+ (GeV) 

FIG. 8: The dependence of the normalized likelihood on the 
decay constant f D +. 
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